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In the study of soft x-ray transitions in solids, there has always been some hope that the results pro- 
vide a direct measure of the density of states. This assumes that (a) matrix element variations over the 
band and (h) final state interactions are small. Both of these assumptions are now known to be incorrect. 
To illustrate the possible strength of these effects, two approximate calculations are presented: the one 
electron oscillator strength of a simple bcc metal as a function of energy; and the strength of the 
Nozieres — DeDominicis singularity at threshold, with phase shifts estimated from an assumed Yukawa 
interaction between conduction electrons and core hole. 
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I. Introduction 

The absorption of a photon can cause an electron to 
change its state. The traditional viewpoint of this transi- 
tion assumes that if one knew the electrons initial i//, 
and final i/// state, then the oscillator strength was 
simply proportional to the square of the momentum 
matrix element. 
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p tf = d^tpipi 



(1.1) 



This simple viewpoint is now known to be incorrect. 
The proper picture is that the optical transition creates 
an electron and hole, and these two excitations interact 
with each other, and each separately with their environ- 
ment. The coulomb scattering between the electron 
and hole is called the exciton effect, named after the 
first known example of Frankel excitons in solids. The 
electron and hole can also emit phonons, suffer elec- 
tron-electron collisions, etc. The sum of all of these 
processes are called final-state interactions. The rate 
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of optical absorption is affected by these subsequent in- 
teractions of the electron and hole. 

Upon reviewing our present understanding of the op- 
tical properties of solids, one finds that some solids are 
well understood while others are not. We classify as un- 
derstood the semiconductors Si and Ge, [ 1] and simple 
metals like aluminum [2,3]. It is noteworthy that in 
these materials the final-state interactions are small: 
they are small in semiconductors because the large 
dielectric constant suppresses the exciton effect. 
Ehrenreich and his co-workers have shown that many 
body effects are small in aluminum. In these solids the 
simple one-electron picture implied in (1.1) seems to 
work quite well. 

We classify wide band gap insulators [4], and metals 
such as copper [5], as solids whose optical properties 
are not well understood. There is not yet good agree- 
ment between theory and experiment in these materi- 
als. The dominant optical transitions in these materials 
create hole states which are heavy, and which are sub- 
ject to significant final state interactions [3]. We 
speculate that final state interactions are important in 
these materials, and explain the difference between 
theory and experiments. 

In the study of x-ray transition in solids, there has al- 
ways been some hope that the results provide a direct 
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measure of the density of states [6]. This presumes 
that (a) matrix element variations over the band are 
small, and (b) final state interactions are small. Both of 
these assumptions are now known to be incorrect 
[7,8,9]. 

The importance of exciton effects in x-ray transitions 
was reported in an earlier reference [10]. The effects 
are large if the hole is highly localized, and if the con- 
duction band of the electron is isotropic. Because the 
coulomb scattering between the electron and hole oc- 
curs in an electron gas, the electron can only scatter 
into states not already occupied. This Fermi-Dirac ex- 
clusion, as well as the exchange interaction among the 
electrons, causes the problem to resemble the Kondo 
effect more than the Wannier-exciton case. The hole 
also has a large interaction with the entire electron 
gas — this leads to a renormalization effect which in- 
hibits the x-ray transition near threshold [11]. Nozieres 
and DeDominicis [ 12] have shown that the exciton and 
renormalization effects combine to give a threshold 
behavior for e 2 (w) 
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where ojt is the absorption threshold frequency, £ ~ 
Ep is a characteristic band energy, and ai is a function 
of the phase shifts di of electrons at the Fermi energy 
scattering from the static potential of the hole. 



«, = 1*-22(2Z+1)(S 



(1.3) 



The term 28ilrr is the exciton part which tends to make 
the threshold singular, while the second term in (1.3) 
arises from renormalization. The angular momentum / 
is that of the conduction electron [13] . If the x-ray hole 
has s-symmetry, the conduction election must have p- 
symmetry (/ = 1). If the hole has p-symmetry, the con- 
duction electron can have either s- or d-symmetry; in 
this case (1.2) has a separate term for each symmetry 
type. 

The x-ray spectra can only be unravelled with a 
knowledge of the phase shifts 8/ — these are the phase 
shifts for simple electron-hole scattering. These have 
been calculated for a free electron gas by Ausman and 
Glick [13]. They find a > 0, and <x x < for / ^ 1 so that 
singularities only occur for the 1 = case (p-state hole). 
These results qualitatively agree with the experimental 
results. 

We have independently calculated the phase shifts, 
and we describe our result in section II. We have con- 
cluded that these phase shifts are qualitatively correct 



but for the wrong reasons. The wrong potential is used 
in the calculation, but it does not seem to make much 
difference in this case. 

II. Phase Shifts 

The present calculations have been performed by as- 
suming that the screened electron-hole interaction has 
the Yukawa form 



V(r) 



exp (— k s r) 



where the Fermi-Thomas screening length is 

k$ = 67Te 2 no/E F 

for an electron gas of density m and Fermi energy E f . 
This is only a crude approximation to the actual poten- 
tial the electron feels when scattering from the hole. 
For example, studies of a point charge impurity in an 
electron gas show that its potential differs from a Yu- 
kawa at long range where Friedel oscillations occur, 
and at short range where its potential is less steep [ 14] . 
In addition, because the hole is part of an atom, there 
are term due to the exchange and orthogonality with the 
atomic wave functions [7,8] . In a pseudo-potential for- 
malism these latter effects contribute a short-range 
repulsive term to the interaction. In spite of these short- 
comings of the Yukawa potential, we believe that it pre- 
dicts phase shifts which are qualitatively correct. The 
reasons for this will be presented below. 

The phase shifts are defined in terms of the eigen- 
values E k and eigen functions ijJk(r) of the electron in 
the region of the potential. 



2m 



+ V(r)-E k )Mr)=0 



(2.1) 



The wave is decomposed into spherical harmonics 

«Mr) = 2 (2l-hl)fP l (k'?)d> l (k,r) (2.2) 

1=0 

In solving for the eigen functions in (2.1), one has the 
choice of specifying the boundary conditions for states 
which are plane wave-like outside of the potential re- 
gion. By choosing standing wave conditions, one is solv- 
ing for a reaction matrix Ki(k,k f ) and the phase shifts 
are denned in terms of the diagonal k = k' component 
[15]. 

tan 8i(k)=-2mkKi(k,k) 

=-2mk l*i'drji{kr)V{r)<ln(k,r) (2.3) 
Jo 
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Similarly, if one chooses outgoing wave conditions, 
then one is solving for a T-matrix ti(k,k') whose 
diagonal components give the phase shift [15] . 

sin 8ie' a i= — 2mkti(k,k) 

= -2mk | r*drji{kr)V(r)<l>i{k,r) (2.4) 

I 
I 

We have chosen to use the reaction-matrix formalism 
(2.3), mostly because it is a real function and this sim- 
plifies computation. 

The phase shifts 8i(k) go to zero as &— » °°. Let us now 
examine their behavior as A ^ 0. From (2.3) we have 
j t (kr) ~ (kr)\ so 



tan8,(£) ~k l + x 



J> 



drV(r)(f),(0,r) (2.5) 



From Levinson's Theorem we know that 



5/(0) = ttMi 



where Mi is the number of bound states of (2.1) with 
angular momentum /. For example, if Mi = we get 
\&i~k'+\ Whereas ifMj=l we get 8 t =7T+ck l+l (where 
c is some constant) for small k. A typical case is shown 
in figure 1, where M =l, and M, = for / ^ 1. So the 
s-wave phase shift comes into tt with a linear slope, 
while all other 8/~ k l+l at small k. 

So in calculating phase shifts, the first thing to de- 
cide is the number of bound states M\. This is deter- 
mined by examining the radial part of Schnidinger's 
equation (2.1) which we put into dimensionless form 
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In (2.7) we have written \ in terms of the density 
parameter r s for an electron gas. Since metallic densi- 
ties vary between 2 < r s < 6, then the range of X values 
in metals is 1.8 < A. < 3.1. For an electron gas of 
metallic density, there is one bound state for s-states, 
and no bound states for any other value of angular 
momentum /. 

These are the predictions of the Yukawa potential. 
We must decide whether these are reasonable conclu- 
sions for the problem of interest. For an actual point 
charge in an electron gas, e.g. , a proton, there is 
probably a bound state. For an impurity in a host metal, 
e.g. , Al atoms in Mg, there is probably not a real bound 
state. This is because the atomic core of the impurity 
cuts off the attractive potential in the region where it is 
the strongest; for example, see Ashcroft's remarks 
about the cancellation influence of the atomic cores 
[8]. 

In the x-ray problem there is certainly an atomic 
core. Yet there is also certainly a bound state in the 
potential. This bound state is the x-ray level itself. That 
is, if an electron scattering from the potential of the 
hole did not think a bound state existed in the potential, 
then it would have no inclination to recombine with the 
hole in the final emission process. Since the emission 
can occur, a bound state does exist which must be 
reflected on the phase shifts. 

The cancellation of the potential in atomic core is 
caused by the necessity of the conduction electron 
wave function to be orthogonal to all of the core states. 
In the x-ray problem, one core electron is absent so that 
the cancellation requirements are less stringent. 

The foregoing discussion shows that any phase shifts 
calculated from a simple Yukawa potential are only 
going to be qualitatively correct. Yet there is some in- 
terest in what this simple model predicts. Our method 
of calculation proceeds by solving directly the scatter- 
ing equation for the reaction matrix. 



K,(k 1 ,k 2 ) = V l (ki,k i ) 



(2.7) 
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(2.8) 



The parameter X determines the strength of the poten- 
tial. Shey and Schwartz [16] have computed the 
number of bound states which exist for each value of \ 
and /. For s-waves (/= 0), they find no bound states 
exist for\<1.68, one bound state for 1.68<\<6.45, two 
bound states for 6.45<\< 14.3, etc. For p-waves (/= 1) 
bound states only exist for A>9.08, and d-wave bound 
states exist for A>21.8. 



where the I th component of the potential V\ is obtained 
from the Fourier transform F(ki-k 2 ) of V(r) 

Vi(k u k 2 )=l ^ dO sin <9P/(cos 0)F(k,-k 2 ) 
£ Jo 



2k, k 2 



Qi({k\ + kl + k*)l2kik 2 ) (2.9) 
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when cos 6 is the angle between ki and k 2 , and Qi is a 
Legendre function. We will abbreviate the argument of 
the Legendre function to write it as Qi(k u A 2 ). Noyes' 
method is used to evaluate (2.8). We find 

Ki(k, k) = Vt(k, k)/(l-A,(k)) (2.10) 

TiViyk, k) Jo p- — k- 
f l (k,p) = V l (k,p) 



+ 



4<m f a 
irF«(A,A) Jo 



p f2 dp'fi(k,p f ) 



ir,(p,t)Fi(*,p') 



-Vi(p,p')Vi(k, k)] (2.12) 



One obtains fi(k,p) from (2.12), perhaps by iterating 
this equation. This is used in (2.11) to obtain A/(/c), and 
thus one has the reaction matrix in (2.10). This is an 
exact result if fi(k,p) is found exactly. The Born Ap- 
proximation result is obtained by setting Ki(k, k)= 
Vi(k, k). 

Note that it is natural to write A/( k) as a power series 
in the interaction strength A= 2/ ksCtsin (2.7) 



A,(*)=_£ ^ m g\ m) (klks) 



(2.13) 



Each successive power of \ corresponds to another 
iteration of the equation (2.12) which determines 
fi{k,p). For example, the first term is 

For / = we can evaluate the integral and express the 
integral as a summation. 



o(l) 

#o 



{k) =mkkji V" i[2p ' sin(W) 



(2.15) 



-sin (2/<p)] 
p=[l + 4/J' 2 /^]-" 2 

(p= tan _1 (2/7A" s ) 
()„(/•. k) = | In (l+4/' 2 /A- 2 ) 

This form is convenient for numerical computation. In 
figure 1 is shown the s-wave phase shift calculated by 
approximating A<> by the first term in (2.13) 

Qo(k,k) 



tan 8o{k) = 



ka B [l-\g[»(k)] 



Also shown in figure 1 are the phase shifts 3\ and 8 2 . 
We also calculated the first correction term gi {1 \k) to 
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Figure 1. The phase shifts 8/ calculated by a Yukawa potential. 
This potential represents a coulomb interaction of unit charge, 
with Fermi-Thomas screening. The phase shift 8 2 is calculated in 
the Born Approximation, while bo and 8\ have corrections for 
multiple scattering. 

the p-wave phase shift, and this correction has been 
added into the calculation of the curve labeled Si. This 
only changes the Born Approximation result by 10 per- 
cent. The 1 = 2 phase shift is the Born Approximation 
result. The multiple scattering terms are small except 
for the s-wave. 

In the x-ray transition, the singular effects occur at 
the Fermi surface, so we are interested in the phase 
shifts evaluated at kj. Figure 2 shows the critical ex- 




FlGURE 2. The variation with electron density r x of the exponents a 
and oli. These depend upon the phase shifts evaluated at the Fermi 
surface. Also shown is the renormalization parameter e and the 
Friedel sum rate result Z. The value of Z is rigorously unity, and 
our deviations from that value indicate the errors in the calculation. 
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ponents a and cx u in (1.3) plotted versus electron 
density /. s . Values of oti for / ^ 2 are essentially equal to 
— 6. Also shown are the values of the Friedel sum Z 
and the Anderson exponent €, where 

1 = 1 J (2/+1)(8,/tt) 



e=2£ (2Z+1) {SiIttY 
i=o 

The Friedel sum rule should rigorously be given by 
z= 1. Our deviation ~20 percent from the exact result 
of z = 1 provides an estimate of the error in the calcula- 
tion. 

Another estimate of the accuracy and convergence 
of the calculational method is obtained by seeing how 
well it estimates the position of bound states. As k-~* 
we get 



im 8o(A) = tan - 



2k/ks 
l-A () (0) 



Abound state is predicted if A >1. We find that 
^(0) = l/2 
4 2 >(0) = ln (4/3) -1/4 

If we just use the first term in the expansion (2.13) so 
A (0)=A./2 then the criteria A o (0) >1 means X>2. This 
is an error of 16 percent from the actual criteria 
X > 1.68. But using two terms gives X/2 + X 2 (ln 
4/3— l/4)> 1 or X> 1.765. So a considerable improve- 
ment in accuracy is obtained by including the second 
term in (2.13). So we conclude that the series (2.13) con- 
verges rapidly for the range of X values of interest. 

III. Matrix Elements 

The change of the matrix element with energy has 
been discussed before and is known to be a large effect 
[7,8,9]. We have calculated the change in matrix ele- 
ment near the first critical point in a bcc solid for an x- 
ray transition from a Is core level. We used a simple 
two band model in each 1/12 of the Brillouin Zone [17] , 
so that the energies and eigenstates are 



JS*=|(8r*+*kHs): 



(g'fc — ^k-G) 2 + ^ 2 



!/2 



|[i(r fc -r k - G ) 2 +P G 



N k - 



The matrix elements were calculated assuming that the 
core was a delta function, and no attempt was made to 
orthogonalize the conduction states to the core. Thus, 
after averaging over polarizations, the matrix elements 
are 

<pl>-=i^v|{^l^*-^*l+(k— G^iar* 

-£*|±2k-(k-c;)JM 

In figure 3 we show the density of states p(E), and also 
the absorption strength 



A(E) 






;.>+8(£-£ A +)+<pj?.>-o(£-£ A 7)} 



In this figure, energy has been normalized to E = 
h 2 G 2 l8m so that & = EIE and v=Vg/Eo. The choice v 
= 0.2 is close to Ham's value for Li (^ = 0.23). Indeed, 
the present calculation was done with Li in mind, since 
the lack of core orthogonalization should not matter 
here. 
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(a) DENSITY OF STATES 



(b) ABSORPTION STRENGTH 



I/=-0.2 




FIGURE 3. The density of states p(e) and absorption strength A{e) 
for a bcc solid such as Li. The energy scale is listed in units 
Eo = h' 1 G 2 iu)l8m, and v—Vno/Eo. The absorption strength A (e) is 
defined as the integral of the x-ray matrix element <p 2 > averaged 
over the energy band. 

The curve for ^ = 0.2 in figure 3(b) has the expected 
shape for Li [18]. Because the wave function is p-like 
at the lower critical point ^=0.8, the transition is al- 
lowed and the absorption strength has the same shape 
as the density of states in figure 3a. The density of 
states structure at if = 1.2 is washed out in A(«f) because 
at the critical point the transition is s-like and largely 
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forbidden. The curve for v = —0.2 has the same density 
of states, but now the upper critical point is p-like and 
A(if) has the same structure as p(W) near 1.2. Neither 
the curve v = Q.2 nor v = — Q.2 has any striking resem- 
blance to the density of states [18]. 



IV. Discussion 

Exciton effects should influence nearly all parts of 
the absorption spectra. The singular behavior at the ab- 
sorption or emission edge is just one prominent feature. 
Another result of final state interactions is that oscilla- 
tor strength is moved from one frequency range to 
another. Often these shifts are small and can be 
neglected. Yet in most cases the true magnitude of 
these effects are unknown, because the relevant 
theoretical calculations are too complicated to do 
realistically. 

The calculation of exciton effects is a formidable task 
which has not been performed properly. Of course 
Wannier excitons have been studied in detail. Some 
model calculations exist near critical point edges — the 
hyperbolic excitons [ 19,20,21] . But the main optical ab- 
sorption strength comes from states throughout the 
Brillouin Zone [5], and all of these states are included 
in the final state interactions. 

So before x-ray absorption and emission measure- 
ments can be used to provide information on the densi- 
ty of states, two sizable corrections need to be made. 
One of these is the exciton effect, and the other is the 
change in the matrix element with energy. 
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